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We study topological insulators characterized by the integer topological invariant Z, in even and 
odd spacial dimensions. These are well understood in case when there are no interactions. We 
extend the earlier work on this subject to construct their topological invariants in terms of their 
Green's functions. In this form, they can be used even if there are interactions. Specializing to one 
and two spacial dimensions, we further show that if two topologically distinct topological insulators 
border each other, the difference of their topological invariants is equal to the difference between 
the number of zero energy boundary excitations and the number of zeroes of the Green's function 
at the boundary. In the absence of interactions Green's functions have no zeroes thus there are 
always edge states at the boundary, as is well known. In the presence of interactions, in principle 
Green's functions could have zeroes. In that case, there could be no edge states at the boundary 
of two topological insulators with different topological invariants. This may provide an alternative 
explanation to the recent results on one dimensional interacting topological insulators. 

PACS numbers: 73.43.-f, 73.20.-i 



I. INTRODUCTION 

The original proposals establishing existence of topo- 
logical insulators with the Z 2 topological invariant 
(in the presence of spin-orbit coupling with time- 
reversal invariance) 1,2 and their subsequent observations 
in experiment 3 ' contributed to a start of a tremendous 
amount of activity on this subject. 

Recently a classification table was worked out for all 
possible noninteracting topological insulators 5,6 . 10 sym- 
metry classes of topological insulators were identified. In 
each spacial dimensionality representative systems of 5 
classes are topologically trivial, while the remaining 5 
classes can be topologically nontrivial. 

The term "topologically nontrivial" in this context im- 
plies the following. First of all, a noninteracting topolog- 
ical insulator, when coupled to an external electromag- 
netic field, is typically (but not always) expected to have 
a response to this field given by a Chern-Simons term in 
even dimensional space or by a 0-tcrm in odd dimensional 
space 5,7 . Second, a topologically nontrivial insulator has 
zero energy excitations at its boundary (or at the bound- 
ary of two topologically distinct insulators). These are 
usually referred to as edge states. To derive all of these 
properties, typically one takes advantage of the fact that 
these systems are noninteracting fermions and thus are 
essentially exactly solvable. 

A natural question is then whether adding interactions 
to topological insulators preserves their observable prop- 
erties, such as the existence and robustness of the edge 
states. 

One particular approach to this problem was developed 
over 20 years ago in a series of papers by G. Volovik and 
collaborators, summarized in his book 8 . This approach 
is based on the following simple idea. Noninteracting 
fcrmionic systems are described by a Hamiltonian which 
is equivalent to a matrix "H a b(k), which is a function of 
momenta k, as explained in Ref. [5]. This matrix satisfies 



certain conditions, which define the symmetry class of the 
topological insulator. For example, for insulators with 
the chiral (sublattice) symmetry, the condition reads 

H(k) = -Etft(k)£, (1) 

where E is a certain unitary matrix. Whether a particu- 
lar system belongs to a non-trivial topological class can 
be investigated by studying the topological properties of 
the function "H(k) which defines a map from the space of 
the momenta k to the space of the matrices %. 

Interacting systems are no longer described by a matrix 
of this kind. However, one can define a Green's function 
(single particle retarded or advanced Green's function) 
for both interacting and noninteracting systems. For 
noninteracting systems, it takes especially simple form 

G(o,k) = p-nck)}- 1 . (2) 

The Green's function computed at zero frequency is sim- 
ply equal to minus inverse of rl and thus inherits all its 
symmetries, implying in case of Eq. (1), that 

G(0,k) = -£tG(0,k)£. (3) 

At this point, we say that an interacting system belongs 
to a particular symmetry class if its Green's function 
G(0,k) satisfies the relations such as Eq. (3). Subse- 
quently, we can study whether a particular interacting 
system is topological by studying the function G(0,k), 
instead of the earlier H(k). 

All of this would stay in the realm of pure math, if 
it were not for a further observation in Ref. [8], on the 
example of symmetry class D systems in 2 dimensions 
(class D is a shorthand for a superconductor without 
any symmetries, such as 2D p + ip superconductor), that 
two topologically distinct systems (at least in two dimen- 
sional space) must have zero energy excitations at their 
boundary. This observation was made by examining the 
Green's function G(fl, k) only, without any reference to 
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whether the system is interacting or nonintcracting, al- 
though as we will see in this paper, an implicit assump- 
tion of the absence of interactions was indeed made. 

The goal of this paper is to construct the framework 
which could in principle generalize these observations 
to arbitrary number of dimensions for all 10 symmetry 
classes. At the basis of this framework lies the construc- 
tion of the topological invariants in terms of the Green's 
functions in arbitrary number of dimensions. One should 
note that simply rewriting the topological invariants in 
terms of G _1 (0, k) instead of "H(k) is not sufficient. To be 
able to relate the topological invariant to the zero energy 
excitations one needs to have it written in space-time, as 
topological characteristics of the map from fi, k space to 
G. 

Furthermore, in this work we are able to eliminate the 
implicit assumption of the absence of interactions. We do 
it by observing that two topologically distinct systems, 
in the presence of interactions, can have zeroes of the 
Green's functions at their boundary instead of the poles. 
The zero of the Green's function is defined as the zero 
of one of its eigenvalues, or equivalently of its determi- 
nant. If this indeed happens, two interacting topological 
insulators can have no edge states at the boundary be- 
tween them and can be continuously deformed from one 
to another without ever closing the bulk gap. Of course, 
this picture does not explain when this scenario happens, 
rather it just says that it in principle can happen. 

In some cases one can give arguments that zero en- 
ergy excitations cannot simply vanish as interaction is 
increased 9 . These arguments, when they are at work, 
lead us to conclude that the poles cannot simply get re- 
placed by the zeroes. On the other hand, it is also possi- 
ble that while the edge states disappear as single particle 
excitations, they remain as collective excitations and con- 
tinue to carry current. In this case, perhaps it is possible 
to reconcile the picture from Ref. [9] and the one from 
this paper. 

At the same time, recent studies of interacting topolog- 
ical insulators in one dimensional space 10-12 show that in 
some cases the edge states do disappear as interactions 
are increased. To make contact with this work, in this pa- 
per we give the details of the Green's function construc- 
tion in one dimensional space, for systems in symmetry 
class AIII or BDI 5 . 

The rest of the paper is organized as follows. In sec. II 
we explain how to generalize the symmetry classes of 
topological insulators to the case when interactions are 
present. In Sec. Ill we discuss the structure of poles 
and zeroes of the Green's functions. In Sec. IV we ex- 
plain why the zeroes of the Green's function can change 
the topological invariants without any zero energy ex- 
citations. In Sec. V we discuss the application of this 
formalism to two dimensions for systems of the integer 
quantum Hall type, which mostly follows the derivation 
from Ref. [8], but takes into account the possible zeroes 
of the Green's functions. This construction can be gen- 
eralized to any symmetry class without chiral symmetry, 



which are classes A, AI, All, D and C, in even number 
of spacial dimensions 5 . Finally, Sec. VI we study chiral 
topological insulators in remaining classes AIII, BDI, CII, 
Dili, CI, in odd number of spacial dimensions. We espe- 
cially concentrate on one dimensional AIII or BDI topo- 
logical insulators. At the end of the paper in Sec. VII, 
we discuss possible directions of further development of 
this formalism. 



II. GREEN'S FUNCTIONS OF INTERACTING 
TOPOLOGICAL INSULATORS 

We begin by reviewing the symmetries of the nonintcr- 
acting topological insulators. Those are characterized by 
matrix Hamiltonians Ti-ij. Here, and in the rest of this 
paper, the indices i, j denote sites of a lattice, spin, and 
other degrees of freedom present in our problem. The 
total number of values i takes is equal to the number of 
degrees of freedom in the problem, we will denote this 
Df. Thus we take for the Hamiltonian of the problem 

H = a\H ij a j , (4) 

where aj, &i are fermionic creation and annihilation op- 
erators. In this formula, as well as in many others below, 
summation over the indices i, j is implied (in some, but 
not all, cases we will write the summation symbol ex- 
plicitly). We also note that sometimes in what follows i 
will simply mean imaginary unit rather than an index; 
we hope the distinction will be clear from context. 

These may satisfy one of three constraints. The first 
is time reversal invariancc 

U\.U*U T = U. (5) 
The second is charge conjugation symmetry 

UyWUc = -H. (6) 

Here Uc and Ut are unitary matrices, satisfying 

U* C U C = ±1, U* T U T - ±1. (7) 

In other words, each of these two symmetries, if it is 
present, can be implemented in two distinct ways, with 
either +1 in Eq. (7) or with — 1. 

In the presence of these two constraints, the Hamilto- 
nian also automatically satisfies the constrain 

S f -HS = -H, (8) 

where S = U^Ut is a unitary matrix. This is called 
chiral, or sublattice, symmetry. The sublattice symmetry 
can also be present in case when neither time reversal 
invariance or charge conjugation symmetry arc present 
(although it cannot be present if only one of these two 
is missing). All together, taking into account Eq. (7), 
one can see that there are exactly 10 separate symmetry 
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classes of Hamiltonians, distinguished by the presence of 
absence of these constraints 5 ' 6 ' 13,14 . 

Note that the Green's functions, defined by 



(9) 



satisfy the following constraints in case when their parent 
Hamiltonians satisfy either of Eqs. (5), (6), and (8). Time 
reversal invariance leads to 



G(fi) = ulG T (n)u T , 

charge conjugation leads to 

G(Q) = -UlG T (-n)U c , 



(10) 



(11) 



(here the symbol G T denotes the transposed Green's 
function) and sublattice symmetry leads to 



G(fl) = -£tG(-ft)E. 



(12) 



We would like to generalize the symmetries given by 
Eqs. (5), (6), and (8) to the case when interactions are 
present, in such a way that the interacting Green's func- 
tion G(Q) still satisfies either of the relations Eqs. (10), 
(11), and (12). In the rest of this section, we show that 
when properly generalized to the presence of interactions, 
the symmetries of the Hamiltonian lead to exactly the 
same constraints on the Green's functions as in the ab- 
sence of interactions 15 . 

We first need to clarify what we mean by the interact- 
ing Hamiltonian and the interacting Green's function. 
The interacting Hamiltonian is constructed in terms of 
the creation and annihilation operators in the same way 
as the noninteracting one, except it is no longer quadratic 
in these operators. Interacting Green's function is then 
constructed in the following way. 

We introduce the retarded Green's function 16 
It can be analytically continued onto the upper half plane 
of complex 0, where, on the purely imaginary axis, it 
becomes the imaginary time ordered Green's function 
Gij{iu>), where O = iu and uj will be taken to be real 
throughout this text. Since the system is supposed to be 
insulating and has a gap in the spectrum, Gij(Vt) can also 
be continued onto the lower half plane through the gap 
where, just below the real axis, it becomes the advanced 
Green's function. For the real values of VL where the 
density of states of the problem is not zero, the Green's 
function has cuts on the real axis through which it can- 
not be analytically continued (if it is, it will spill into 
another Riemannian sheet). In some cases we will be in- 
terested in topological insulators with boundary modes 
which will manifest themselves in excitations at zero and 
close to zero. In this case, to define the Green's function 
analytic everywhere we can think of a finite system so 
that the singularities of G(f2) are confined to a number 
of poles. We will elaborate on the analytic properties of 
the Green's function in the next section. 

In addition to the Green's functions defined here, it 
is often customary to introduce time-ordered, anti-time 



ordered and Keldysh Green's functions 17 . To distinguish 
them from the Green's functions defined in the previ- 
ous paragraph, we will always specify which function we 
deal with whenever we work with one of those functions, 
retaining the name "Green's function" for the analytic 
retarded-advanced-imaginary-time-ordered functions. 

Now let us generalize the symmetries of the Hamilto- 
nian to the case when interactions are present. We define 
unitary operators Ut, Uc and E, which act on creation 
and annihilation operators a] and a, (where i labels the 
sites of a lattice as well as other degrees of freedom if 
present) in the following way 



j 



UTij&j, U±a}U T = 



UtfT„-X^ & \ UT 



j 



t 



(13) 



U^afic = ]T tf Cij a], tf c a\Uc = £ a 3 U c ^ (14) 



and 



t^a t t = ]T a]E]. i; Stats = £ s tf&j .. (i 5) 



We now declare that a many body Hamiltonian H is time 
reversal invariant if 



(j\.H*U T = H, 
is particle-hole symmetric if 

ulHUc = H, 

and is chiral if 

£*H*£ = H. 



(16) 



(17) 



(18) 



Here H* is a complex conjugate many body Hamiltonian. 
The creation and annihilation operators are not affected 
by the operation of complex conjugation. 

Now it is straightforward to check that if the Hamilto- 
nian is noninteracting, that is if 



H = ^ a\%ija,j, 



(19) 



then Eqs. (16) and (13) are equivalent to Eq. (5), 
Eqs. (17) and (14) arc equivalent to Eq. (6), while 
Eqs. (18) and (15) are equivalent to Eq. (8). For ex- 
ample, 



Uj,H*Ut = V\a\U*{afi T = a\u\, ki U* i UT i iai = 



aiHij&j = H. 



(20) 



Similarly, 



U^HUc = Ul&lUijajUc = a k U C kt'H lj ul, :jk al 



(21) 
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(this assumes, without loss of generality, that trH = 0), 
and finally 

-alTtljUjiSikak = a}Uija,j = H. (22) 

The generalizations given here are not the only way to 
generalize these symmetries to the case when interactions 
are present. One alternative way to generalize the chiral 
symmetry is given in appendix A. However, they seem 
to be most natural from the point of view of realistic 
applications. 

Next, one can check that even in the presence of in- 
teractions Eqs. (10), (11) and (12) still hold. This is 
a straightforward yet tedious calculation which can be 
skipped on the first reading. 

We work with the time ordered and anti-time ordered 
Green's functions 

Gjjit) = -i(0\Tai(t)a](p)\0), 

G\ 3 {t) = -z<0|T-a,(t) at(0)|0>. (23) 

Here the symbols T and T denote time and anti-time 
ordering of subsequent time dependent operators. We 
denote the corresponding Green's functions G* and G* 
(to avoid confusion with the notation G T which implies 
transposed Green's function Gfj = Gji). Given these 
time ordered Green's functions, we can construct re- 
garded and advanced Green's functions according to 

Gf 3 (t) - 6(t) (G^t) - G^t)) , 

Gfct) = 0(-t)(G^.(t)-4.(t)). (24) 

Fourier transform of those gives us cither G R (ft) with 
complex ft in the upper half plane or G A (Cl) with com- 
plex fl in the lower half plane, which are known to be 
analytic continuation of each other (at least in the case of 
the gapful spectrum), and define the same analytic func- 
tion G(f2) in the entire complex plane. To prove that thus 
defined G(f2) satisfies Eqs. (10), (11) or (12), depending 
on the symmetries present, it is sufficient to prove that 
the time ordered and anti-time ordered Green's function 
satisfy one of the three relations, 

G\t) = u\,G tT (t)U T , (25) 

G\t) = -tf c & T {-t)V c , (26) 

and 

G\t) = -£tG*(-t)E (27) 

(same for the anti-time ordered Green's function). 

Start with the time reversal invariance. The presence 
of the symmetry Eq. (16) implies that for every many- 
body state \n), there exists a state Uj,\n*) which has 



exactly the same energy as the state \n) (here \n*) de- 
notes the complex conjugate state). Indeed, assuming 
that H\n) = E n \n) gives 

M f T \n*) = ulu T HU f T \n*) = U^H* \n*) = E n U f T \n*) . 

(28) 

Without loss of generality we assume that the ground 
state is time reversal invariant, that is 

III |0*> = |0) , |0*> = (It |0) . (29) 

If it were not time reversal invariant, we could still take 
expressions such as |0)±C4|0*) as a ground state, which 
arc invariant. 

Now assume that t > 0. Then the following sequence 
of transformations can be carried out 

iG\ j {t) = {Q\e i " t a i e- i " t a)\Q) = 

(0| U^e^UT^U^e-^UTa^U^UT |0) = 

(0*| e^^UTCLiU^e-^^UTa]^ |0*) = 

U Tik (O'le^ake-^&l^Uij = 

4 ifc (0\a l e- i6t aie i6t \0)U lj = 

tf Tik iG tT kl (t)UTir (30) 

To go from the third line to the fourth line of Eq. (30) 
we took advantage of Eq. (29). To go from the third to 
the fourth line, we used that for general states |n) and 
\m) and for generic operator A, the following is obviously 
true 

(n\A\m) = (m*\A T \n*}. (31) 

We could also prove a similar relation for t < 0, but this 
is not really necessary since we can just use the retarded 
Green's functions to construct G(Q). At the same time, 
we can easily prove the same for the anti-time ordered 
function G*. Thus we have proved Eq. (25) and, by ex- 
tension, Eq. (10). 

We now move on to charge conjugation symmetry. Just 
as before, each excited state is doubly degenerate, and 
the ground state satisfies 

U c |0) = |0) . (32) 

We now do the following sequence of transformations for 
t > 

iG%{t) = <0| e^Oie-'^ot |0> = 

(0| Ule^UchUle-^UcaXUcUl |0) = 

U Cj i {0\e iAt ale- i6 %\0)tf Cki = 

-U Cj i (0\Ta ie im aie- iHt \0)Ul ki = 

-U C fliG\ k {-t)Xjl kl (33) 

From this Eq. (26) easily follows. In addition, same re- 
lations hold for negative time t < as well as for the 
anti-time ordered Green's functions G*. Thus we have 
proved Eq. (11). 
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Finally, for the chiral (sublattice) symmetry, we repeat 
the same process. We take the ground state to be invari- 
ant under the chiral transformation, which leads to 



S t |0*) = |0), |0*)=E|0). 
We find, for t > 0, 

iG t ij (t) = (0\e iAt a i e- iAt a]\0) = 
(0| ^e^^ta^e-^^taWt |0) 



(34) 



at |0*) Ey = 
Ej fe (0|ate-^a fe e^|0) £y = 
-Ej fc (0\Te- iHt a k e iHt a}\0)X lj = 



(35) 



Eq. (27) follows immediately, and, together with the ex- 
tension of this relation for t < and for G*, Eq. (12) 
follows as well. 

Thus we have proven that Eqs. (10), (11), and (12) 
hold even in the presence of interactions. 

Finally, note that regardless whether or not the system 
is interacting, its Green's function satisfies 

Gy(fi) = G*.(fi*). (36) 

Here is the Hermitian conjugate of the Green's func- 
tion understood as a matrix. 



III. ANALYTIC PROPERTIES OF THE 
GREEN'S FUNCTIONS 

In this section, we discuss the analytic properties of 
the Green's functions, important for subsequent applica- 
tions. The main result of this section is that the deter- 
minant of the Green's function can always be written in 
the following form 



detG 



II^l(fi-*n) 



(37) 



Here r s are real numbers representing positions of the 
zeroes of the determinant, e„ are the excitation energies 
which are also obviously real and which represent the 
poles of the determinant, Dh is the combined dimension 
of the Hilbert spaces of the system with one extra par- 
ticle than the ground state and one less particle than 
the ground state, and Df is the number of degrees of 
freedom, or number of distinct creation or annihilation 
operators as discussed earlier. Generally 



D f < D h 



(38) 



where the equality is achieved, in particular, if the system 
is not interacting. This can be seen if one notes that in 
the presence of interactions consists of all the states 
obtained from the ground state by adding or removing 



a particle (which can be done in Nf distinct ways) plus 
states which differ from those by particle-hole excitations. 

If the system is not interacting, = Df and the 
determinant of its Green's function is particularly easy 
to calculate starting from G = [Cl — 'H]~, to give 



det G=H 



n 



(39) 



which is consistent with Eq. (37). 

In the presence of interactions, the determinant of the 
Green's functions, in addition to its poles, has zeroes as 
follows from Eq. (37). Unlike the poles which have clear 
physical meaning as excitation energies, the zeroes are 
more subtle and, as we will see below, their positions re- 
flect the matrix elements of the creation and annihilation 
operators. However, it should be clear that as interac- 
tions of the system are switched off, the zeroes move to 
coincide with the positions of some of the poles, resulting 
in disappearance of one zero and one pole at a time from 
the determinant, until in the absence of interactions it 
takes the form Eq. (39). Conversely, if the interactions 
are turned on, the zeroes and poles appear in pairs. 

In the rest of this section we present the proof of 
Eq. (37). This can be skipped on the first reading. 
We begin by writing down the spectral decomposition 
(sometimes referred to as Lehmann decomposition) of the 
Green's function 



(40) 

Here \n) are the states which have one more particle com- 
pared to the ground state, and 



,+ - 



= E n (N+l)-E o (N)>0, 



(41) 



where E (N) is the ground state of the system of N par- 
ticles (we take the chemical potential as included in the 
Hamiltonian, thus dE (N)/dN = 0), and E n (N + 1) is 
the excited state of a system with one more particle. Sim- 
ilarly, |m) have one less particle than the ground state, 
and 



u~ =E (N)-E m (N-l)<0. 



(42) 



It is convenient for the moment to work with combined 
notations, where all states are denoted by \n) regardless 
or not whether they have one more particle or one less 
particle than the ground state and all the energy is de- 
noted by e n which can be either positive or negative and 
coincide with cither uj, 



or uj m . Then we can write 



(43) 



where U n j — (n\ Oj |0) if \n) has one more particle com- 
pared to the ground state and U n j = (0| a'- \n) if \n) has 
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one less particle. Note that 



E^ 



(44) 



In the matrix U n j , n goes over exactly Dh values while j 
goes over Df values, where, as we recall, Df < Dh- 

Now we are in the position to study the determinant 
of the Green's function detG, to derive Eq. (37). 

If the system is non-interacting, then it has as many 
excited states as the number of its degrees of freedom 
or Df — D^. Then U n j is a unitary matrix which can 
be diagonalized. It follows that the eigenvalues of the 
Green's functions are 



(45) 



so that we immediately reproduce Eq. (39). 

If, however, the system is interacting, then Dh can 
be larger than Df and U n j is a rectangular matrix. In 
that case, we write down the formal expression for the 
determinant of the Green's function, starting from the 
definition Eq. (43) 



detG =E e n u ^ UnMi) (-i) p(ct) 



(46) 



Here as usual cr(i) is a permutation of the numbers 
1,2,..., Nf, and P denotes its parity. 

It is now easy to see that if any two n coincide, for 
example if n\ = ri2, then the summation over permuta- 
tions gives zero. This means that the determinant of the 
Green's function has simple poles at each point e„ (and 
does not have any higher order poles). This justifies the 
denominator of Eq. (37). 

To explain the origin of the numerator, we first note 
that the determinant of the Green's function is clearly 
a rational function of f2, thus the numerator must be a 
polynomial. Moreover, it is well known 16 that at very 
large O the Green's function asymptotically approaches 
a very simple form Gij ~ 5ij/Cl. This means that its de- 
terminant must go as 1/ il D f . This implies that the poly- 
nomial in the numerator should be of the order Dh — Df, 
with the highest order term Sl Dh ~ D f having the coeffi- 
cient equal to unity. 

Thus the numerator of the determinant must take pre- 
cisely the form shown in Eq. (37). The only thing which 
we have not yet proven is that all r s are real. 

To prove that, let us prove that the eigenvalues of 
the Green's function cannot have zeroes (cannot vanish) 
when f2 is in the upper half plane or in the lower half 
plane. In this proof we follow the arguments given in 
Refs. [16,18]. Let us examine the eigenvalue equation 



^^k^) (0) = A ^ m w (Q) . (47) 



jn 



Then for Q — E where E is a real number which makes 
the Green's function hermitian 



w = E 



E-e n 



E 



Z, 



E-e' 



(48) 



where the notation 



> 



(49) 



was introduced. In addition, if ft is taken to infinity, 
A P (S7) ~ as we have already discussed. 

Now let us take ft around the contour shown on Fig. 1. 
The horizontal part of the contour goes infinitesimally 
above the real axis. On this contour, it follows from 
Eqs. (48) and (49) that the imaginary part of an eigen- 
value A M (f2) is always zero or negative (negative if Q is 
close to one of e n ). At the same time, when f2 goes over 
a large semicircle, A M (S1) = and its imaginary part 
is also negative. Thus ImA M (fi) < for all il belonging 
to the contour shown on Fig. 1. 

However, there exist a well known theorem in complex 
analysis which states that if an analytic function X(il) 
has poles or zeroes in a certain domain, then if one plots 
A(f2) as encircles this domain counterclockwise, the 
number of times A(f2) encircles zero in the counterclock- 
wise direction is equal to the number of zeroes minus the 
number of poles of this analytic function inside this do- 
main. We see that as CI follows the contour shown on 
Fig. 1, A /i (Q) does not encircle zero at all. Since A p (f2) 
cannot have any poles in the upper half plane, it follows 
that it cannot have any zeroes in the upper half plane. 

Similarly, A M (f2) cannot have any zeroes also in the 
lower half plane. This means they can have zeroes only 
on the real axis. 

If the eigenvalues of a matrix have zeroes on the real 
axis only, this implies that the determinant of the matrix 
can have zeroes on the real axis only as well. Then the 
fact that all r s in Eq. (37) are real follows. 

This concludes our proof of Eq. (37). 



IV. MODE COUNTING VIA THE GREEN'S 
FUNCTIONS 

At the heart of the approach to the topological insula- 
tors developed by Volovik 8 lies the following topological 
invariant 



Aq = Tr 



duj 



iri 



- Kd^G. 



(50) 



Here, first of all, a; is a real variable such that = iuj. 



Second, is a matrix inverse to Gij 



(51) 
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the spectrum). If the system is noninteracting, the only 
way it can increase the number of positive energy levels 
is by moving an energy level from negative to positive 
values (and thus closing the gap somewhere in the mid- 
dle of the process where this energy level crosses zero). 
If the system is interacting, it can change the number 
of positive energy levels by moving a zero from positive 
to negative values (where it eventually may "fuse" with 
a pole and remove it from the spectrum). This process 
happens without ever closing the gap. A possible sce- 
nario outlining the motion of zeroes and poles is shown 
on Fig. 2. 




FIG. 1: The contour Q is taken along in order to prove that 
the eigenvalues of the Green's functions cannot vanish in the 
upper half plane. 



The notation G _1 , which would seem to be natural to 
use instead of K in this instance, is going to be reserved 
for a different quantity which we will define later. 

The existence of this topological invariant reflects the 
fact that TTi(GL(Df,C)) = Z. It is straightforward to see 
that if the system is perturbed and its Green's function 
changes by small amount as G — » G + SG, this quantity 
remains unchanged. 

To see the meaning of this quantity, we express it in 
terms of the eigenvalues of the Green's function 



-( 



duj d In (A M ) 



7TZ 



f 

J — c 



duj d In dct G 



duj 



(52) 



It is easy to see by direct substitution of Eq. (37) that 
the invariant is equal to 



(53) 



that is the difference between the total number of pos- 
itive and negative poles, minus the difference between 
the total number of positive and negative zeroes of the 
determinant of the Green's function. 

For noninteracting systems, Ni is simply the difference 
between the number of positive and negative eigenvalues 
of the system 8 , as noninteracting systems have no zeroes. 
For interacting systems, N\ turns out to be more compli- 
cated. Notice, however, that N\ can change only if a pole 
crosses zero (it increases by one if the pole moves from 
negative to positive values), or if a zero crosses zero (it 
increases by one if a zero moves from positive to negative 
values). 

This immediately leads to the following conclusion. 
Suppose we have an insulating system (it has a gap in 



FIG. 2: A possible scenario outlining the change in the po- 
sition of the poles (solid lines) and zeroes (dashed lines) as 
a function of some parameter A, as the interaction strength 
is varied, (a) an energy level crossing zero as a function of 
A for a noninteracing system, (b) Two zeroes and and two 
nearby poles appear in the determinant of the Green's func- 
tion Eq. (37) as the interactions are increased (c) zeroes move 
closer to the positions of the original poles of the noninter- 
acting problem (d) zeroes "annihilate" the poles resulting in 
a system with exactly the same change of the topological in- 
variant, but without any energy levels crossing zero energy. 

This observation lies at the core of the possible transi- 
tions between two topologically distinct insulating states 
without closing the gap, in the presence of interactions. 
We will further elaborate on it later in this paper. 



V. TOPOLOGICAL INSULATORS IN TWO 
DIMENSIONAL SPACE 

In two spacial dimensions there are three types of topo- 
logical insulators characterized by the topological invari- 
ant Z. These are class A (an example is integer quantum 
Hall effect), class D (an example is a p + ip supercon- 
ductor) and class C (an example is a 2D superconductor 
which is not time reversal invariant but with spin rotation 
invariance) 5,13 . More generally, representative systems of 
a total of 5 symmetry classes can possess a Z invariant 
in even spacial dimensions. 

The topological invariant which characterizes these 
systems in the context of integer quantum Hall effect 
in the absence of interactions is referred to as TKNN 
invariant 19 . It is usually written as a flux of the Berry 
connection through the Brillouin zone. However it is also 
possible to write it in terms of the Green's functions, as 
was first observed in Ref. [20]. The advantage of the ex- 
pression in terms of Green's functions is that it is then 
straightforward to generalize it to the case when interac- 
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tions are present. 

To write down this invariant, we first introduce the 
Green's function in the momentum space, G a (,(Q, k) (this 
assumes that our system is translationally invariant). 
Here the indices a, b label the remaining degrees of free- 
dom once momentum is introduced (these can be the ba- 
sis of the lattice, spin, and other degrees of freedom) . Wc 
also introduce the notations fc = u, where f2 = iiv and w 
is real, in addition to the momentum vector components 
k\ and fc 2 . 

Now the TKNN invariant takes the form 20 

N 3 = (54) 

^ptr l°° J ^G-'d^GG-'d^GG-^G, 

where summation over a, ft, 7 taking values 0,1, and 2 is 
implied, and tr denotes the trace over the matrix indices 
of G. G -1 is a matrix inverse of G, defined as in 

Y,G-£{Sl,k)G bc {Sl,k) = 5 ac . (55) 

b 

Following Ref. [8], we denote this invariant by N 3 , to 
emphasize its similarity with Eq. (50). The value of this 
invariant is an integer number which, in case of integer 
quantum Hall effect, is proportional to the Hall conduc- 
tivity a xy . 

However, we can now imagine interactions to our sys- 
tem. The presence of interactions does not change the 
fact that N 3 is a topological invariant. Indeed, its ex- 
istence reflects the homotopy group ir 3 (GL(M, C)) = Z 
for a matrix of a size M by M and is totally unrelated 
to whether the Green's function describes interacting or 
nonintcracting system. 

However, a question still remains whether the existence 
of this topological invariant in a system with interactions 
would have any consequences for its physical properties. 
In Ref. [8] it was shown that in the absence of interac- 
tions, if there are two domains separated by a straight 
one dimensional domain wall where N 3 takes on different 
values on either side of the domain wall, this implies that 
there are zero energy excitations on the domain wall, that 
is, the edge states. 

Let us see what remains of this phenomenon in the 
presence of interactions. We will see that if two domains 
with two different values of the invariant share a domain 
wall which separates them, then even in the presence of 
interactions one can have zero energy excitations on the 
domain wall. However, the new effect which appears in 
the presence of interactions is that instead of a zero en- 
ergy excitation, the Green's function can also develop 
a zero at zero energy on the domain wall. This would 
still lead to the change in topological invariant from the 
first to the second domain. Thus in the presence of in- 
teractions, sometimes it is possible to have to adjacent 
distinct topological insulators without any edge states in 
between. 



To derive this result, we closely follow the discussion 
in Ref. [8]. Let us imagine that we have two distinct 
topological insulators in the domains x > and x < 
of the two dimensional plane (assuming the two dimen- 
sional cartesian coordinates x and y). x = is the po- 
sition of the domain wall. Suppose the parameters of 
the Hamiltonian smoothly change as a function of x, so 
that asymptotically at large positive and large negative x 
the Hamiltonian looks like that of a topological insulator 
with the invariant and N 3 . 

The Green's function of this topological insulator is 
no longer just the function of the momenta, as the mo- 
mentum k x is no longer a good quantum number. How- 
ever, the translational invariance in the y-direction is still 
present. Generally, the Green's function, a matrix in the 
position space, can be written as G a f,(Q; x, x'; y—y')- One 
can perform a Fourier transform of the Green's function 
in the y-direction, along the domain wall, and a Wigner 
transform 17,21 of the Green's function in the x-direction 
perpendicular to the domain wall, to arrive at the Green's 
function G a {,(f2, k, x). 

With the help of the Wigner transformed Green's func- 
tion, we define a local ^-dependent topological invariant 
(fco = oj, k\ = k x , and k^ — k y as before and a, ft, 7 are 
summed over the values 0, 1, and 2) 

N 3 (x) = (56) 

duj / f-^G- 1 ^ GG^d^GG^dk G. 
J (2tt) 2 k " kf> 

Here G _1 (r2,k, x) is defined as satisfying the equation 
Y,G~ b 1 (n, 1 k,x)G bc (Cl,k,x) = 5 ac . (57) 

b 

Here, and below, we use the notation tr to denote the 
trace over matrix indices which do not include coordi- 
nates or momenta, and the notation Tr whenever the 
trace includes the summation over coordinates and mo- 
menta, such as in Eq. (58). 

Crucially, G _1 is no longer the inverse of the Green's 
function and is distinct from the Wigner transform of the 
matrix K defined earlier in Eq. (51). K at the same time 
is the true inverse of the Green's function. Far away from 
the domain wall where the Hamiltonian and the Green's 
function stop depending on x explicitly, G _1 and the 
Wigner transform of K coincide. 

Since the topological invariant can take integer values 
only, as x varies we expect that N 3 charges abruptly from 
N 3 to at some value of x, say x = 0. 

Using a technique developed in Ref. [8], we can show 
that the difference iVjp — N 3 is related the topological 
invariant N\ defined in Eq. (50). To do that, consider N\ 
calculated at a certain value of momentum k y = A > 

/T 
— K{n,ky)d bJ G{n,ky), (58) 
m 

where the symbol Tr denotes summation over all the in- 
dices not represented by p y and the second line of this 
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equation is added for later convenience. Construct the 
difference 



^(tfl|fc„=A- ^iI^-a) 



(59) 



where A is some very large positive momentum. For the 
nonintcracting system, this gives the excess of positive 
energy excitations at k y = A over those at k y = —A, 
as follows from Eq. (53). If this number is nonzero, this 
means that there is a branch of excitations which crossed 
zero as a function of k y . In other words, there are zero 
energy excitations (edge states) at the boundary between 
two topological insulators. 

For interacting system, this is a more subtle quantity 
as we have seen in the previous section. It can be nonzero 
even if there are no edge states, by moving a zero of the 
Green's function across zero energy as k y changes. The 
formalism we are developing cannot tell which of the two 
scenarios happen in general, but nevertheless it can relate 
this number to the difference — N^. 

To derive this, we need a relation between K and G _1 . 
The exact relationship between those is very complicated. 
If however the ^-dependence of the Green's function is 
weak (the Hamiltonian varies slowly as x is taken across 
the domain wall) the relationship between them can be 
derived by the so-called Moyal product expansion or a 
gradient expansion . This is defined in the following 
way. Suppose we have two functions A(x, y) and B(x, y), 
and a third function C(x, y) which is a convolution of the 
first two, 



C(x,z] 



-I 



dyA(x 1 y)B{y 1 z). 



(60) 



Then the Wigncr transform of these functions is approx- 
imately related to the Wigner transform of the first two 
by the following expansion 



C{k,x) = A(k,x)B(k,x) 
1 



(61) 



— (d x A(k, x) d k B(k, x) - d k A(k, x) d x B(k, x)) + 

valid when x dependence of these functions is weak. 

Note that the trace of C can be calculated in a straight- 
forward way 

JdxC(x,x)=J ^C(k,x) = J ^A(k,x)B(k,x). 

The inverse Green's function K's convolution with the 
Green's function G is a delta-function. This allows us to 
calculate K approximately via G _1 with the help of the 
gradient expansion 



K = G -1 + (G~ 1 d x GG~ 1 Gk x GG~ 1 - 

G- 1 d k hc- 1 d x GG- 1 ) + .... 



(63) 



We use this to calculate Eq. (58) using the gradient 
expansion. We substitute Eq. (63) into Eq. (58) (there is 



no need to use Moyal product expansion in the product 
of K and d u G, due to Eq. (62)). We find 

N 1 {k y ) = tx J dx ^^(G- 1 d x GG- 1 d k:c GG- 1 d L0 G- 
G- 1 d kx GG- 1 d x GG- 1 d 0J G) = 

E^tr / d ^^G- 1 d a GG- 1 d p GG- 1 d 1 G. (64) 



a/3 7 



Here a, j3, 7 go over x, k x ,w. 

The next step relies on the construction of a vector 

n a = e aMS ti G- 1 df j GG- 1 d 1 GG- 1 d s G, (65) 

where a, /3, 7, 5 now span the four-dimensional space of 
x, k x , k y , and lj. This vector satisfies 



d a n a = 0, 



(66) 



as can be checked by direct differentiation. The relevant 
difference of the invariant Eq. (59) can be thought of as 
the flux of this vector through a three dimensional surface 
formed by the planes p y = A and p y = —A, surrounding 
the singular point lj = 0, x 
tuned to appropriate values. 

That same flux can be written as a flux through the 
surfaces x = L and x = — L. Finally, this can be brought 
to the form 



0, while k x and k y are 



1 



(i\Ti(A) - JVi(-A)) = N 3 (L) - N 3 (-L). 



(67) 



This formula relates what happens at the boundary to 
the difference of topological invariants. In case of non- 
intcracting system, it proves that for two topologically 
distinct systems bordering each other, there are always 
boundary zero energy excitations (chiral edge states) . In 
the presence of interactions, it is possible for the left hand 
side of this formula to be nonzero without having any 
zero energy excitations at the boundary as we saw be- 
fore. Whether this happens in practice requires further 
study of the zeroes of the Green's functions in interacting 
topological insulators. 



VI. TOPOLOGICAL INSULATORS WITH 
CHIRAL SYMMETRY 

We now proceed to generalize the previous discussion 
to the topological insulators with chiral or sublattice 
symmetry. These include the insulators in classes AIII 
(sublattice symmetry only), BDI (sublattice symmetry 
and time reversal invariance), CII (sublattice symmetry 
and spin orbit coupling), Dili (time reversal invariant 
superconductors with spin orbit coupling such as phase 
B of He 3 ) and CI (spin singlet time reversal invariant 
superconductor), all of which can possess a Z topologi- 
cal invariant in odd spacial dimensions only, unlike the 
remaining insulators including those discussed in the pre- 
vious chapters which could possess a Z topological invari- 
ants only in even spacial dimensions. Expressions for the 
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topological invariants such as Eq. (54) work only in the 
absence of chiral symmetry (they are designed to work in 
even number of spacial dimensions). In what follows we 
propose equivalent expressions for chiral topological insu- 
lators and proceed to relate those to zero energy bound- 
ary excitations and zeroes of the Green's function. 

The simplest of these insulators are the ones in class 
AIII, since they do not have any symmetries besides the 
sublattice symmetry. Let us first review the noninteract- 
ing topological insulators of this type. 



A. Noninteracting chiral topological insulators 

Consider a noninteracting topological insulator in class 
AIII. A noninteracting fermionic problem is described by 
its Hamiltonian 



(68) 



where the indices i, j refer to points on a lattice, spin, 
and/or species of fermions in the problem. Such a system 
belongs to class AIII if there exist a unitary matrix E,j , 



^ij^jk — $ik, 



such that 



E 

jk 



-Hi 



(69) 



(70) 



This Hamiltonian is said to be insulating if it has a single 
particle gap. 

It is straightforward to see that . EjjEjfc = Sik as 
a consequence of Eq. (70), thus the eigenvalues of E are 
cither +1 or — 1. 

The existence of this "chiral" or "sublattice" symme- 
try given by Eq. (70) immediately leads to a number of 
consequences. For any single particle state fa of energy 
e„, satisfying 



Hijipj = e„V»7 



(71) 



there exists a state ^ijipj with energy —e n . Among 
the states there may be zero energy single particle states 
satisfying 



(o) 



0. 



(72) 



Those are also eigenstates of Ey with the eigenvalues ei- 
ther +1 or — 1. The zero energy eigenstates whose eigen- 
value of Ejj is +1 are called the right zero modes, and 
those with the eigenvalue —1 the left zero modes, 



e4 0) 



(73) 



change HWij is deformed smoothly). In fact, that differ- 
ence is called the index of the operator V a p from Eq. (75) 
and the fact that it is a topological invariant is referred 
to as the "index theorem". The existence of the index 
theorem for the chiral topological insulators, as well as its 
absence in non-chiral topological insulators, in particular 
in p + ip superconductors, was emphasized in Ref. 22. 

It is convenient to work in the basis of the eigenstates 
of E. In this basis, E looks like 



5 a f) 






Sa/3 



(74) 



We often refer to the sites of the lattice which arc +1 
eigenstates of E as even and —1 eigenstates as odd sub- 
lattice. The Hamiltonian takes form 



H = 





vl 







(75) 



Here a and f3 label sites of each of the sublattices. It is 
straightforward to check that Eq. (75) is consistent with, 
and follows from, Eq. (70). 

For a translationally invariant system, we can intro- 
duce momentum (crystal momentum if on the lattice). 
In this basis, the Hamiltonian takes the form 



n 







V ab (k) 




(76) 



Here a and b label remaining states within the lattice 
basis after the momentum is introduced. 

Such a Hamiltonian is known to be characterized by a 
topological invariant in any odd dimensional space. The 
invariant is constructed the following way. V a fj is a ma- 
trix function of momentum k. We can construct an ex- 
pression 



f d d k f[ v- 1 

•* n=l 



d n V. 



(77) 



Here d a — dk a and appropriate summation is implied. 
The integral goes over all the entire (i-dimensional space, 
or, in case if the problem is defined on a lattice, over the 
Brillouin zone. This expression takes on values propor- 
tional to an integer, which, mathematically speaking, re- 
flects the fact that n d (GL(M, C)) = Z for odd d. This in- 
variant is the odd-dimensional counterpart of the TKNN 
invariant, which exists in the symmetry classes without 
chiral symmetry only, and is nonzero only in even d. 

It is especially simple to see that this is an integer if 
the matrix V is a one by one matrix, or reduces to a 
number z(k), and if d = 1. Then 



h 



dk ( 1 dz 
47r V z dk 



1 dz* 
z* dk 



(78) 



The difference of the number of the right and left zero 
modes is a topological invariant (this difference cannot 



This measures the number of times z(k) winds around 
the origin as k goes over the Brillouin zone. 
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One famous example of the topologically nontrivial in- 
sulator in one dimension with the symmetry Eq. (70) 
are the Su-Schricffcr-Heeger 23 solitons. Indeed, these au- 
thors considered a problem of a particle hopping on a ID 
lattice with the Hamiltonian 

H = J2(t 4i«2i+i + t' a 2i+1 a 2i+2 ) + h.c. (79) 

i 

It is straightforward to see that the matrix "H in this case 
takes the form 



(80) 



We sec that t + t'e lk winds around the origin if t' > t and 
does not if t' < t. It immediately follows that if for i > 0, 
t > t', and for i < 0, t < t' , at the boundary between 
t > t and t < t' insulators there is going to be a zero 
mode, the Su-Schricffcr-Heeger soliton (this can of course 
be established also in a direct fashion by constructing the 
solutions to the Schrodinger equation). 

Technically speaking the Hamiltonian Eq. (79) pos- 
sesses not only chiral symmetry Eq. (70) but also time- 
reversal invariance, so it belongs not to the class AIII, but 
rather to the class BDI. However, in ID the distinction 
between these two classes is not significant (although it 
becomes more significant in spaces of higher dimension- 
alities). 

More generally, any problem where a particle hops on 
a bipartite lattice possesses the symmetry Eq. (70) and is 
said to belong to class AIII (in the absence of any other 
symmetries), with the matrix X acting by +1 on even 
and by —1 on odd sites, as long as a particle can only 
hope from an odd site to an even site or vice versa. 



B. Interacting chiral topological insulators 

We would like to see if adding interactions to the 
Hamiltonian (68) preserves its topological properties, es- 
pecially the edge states. To do that, we need to generalize 
Eq. (70) to interacting systems. We do it in the way dis- 
cussed earlier in Sec. II, in particular in Eqs. (15), (18) 
and (12). 

As a result of the symmetry defined by Eq. (18) for any 
state \n) there exists a state fit \n*) which has exactly the 
same energy as the state |n), 



H\n)=E \n) , \n*) = E±^ \n*) . 



(81) 



Notice also that if the state \n) has N particles, then its 
conjugate state £t \n*) has Df — N particles. To see that, 
we write 



D f 



D f 



^2&\ait^\n*) = SfE^atoiSt \n*) = 

i=l i=l 

D f / D S 

St£a4K) = St \d s - J^al&i |K). (82) 



The lowest energy state of the system obeying this kind 
of symmetry has exactly N = -D//2 particles, that is, the 
system is at half filling. Indeed, suppose the number of 
particles in the lowest energy state is not Df/2. That 
means that the ground state energy at particle number 
fixed at N, E (N), has a minimum at N ^ -D//2. At 
the same time, it satisfies 



E (N) = E {Df-N) 



(83) 



as a consequence of Eqs. (81) and (82). Therefore, 
this function also has exactly the same minimum at 
N± = Df — N particles. This implies that the compress- 
ibility of the system, d 2 E /dN 2 , is not positive definite. 
From now one we restrict our attention to systems with 
positive compressibility, which as we just saw is equiva- 
lent to systems at half filling. 

An interacting system invariant under Eq. (18) pos- 
sesses a topological invariant, a generalization of Eq. (77) 
for the case when there are interactions and reducing to 
it in case where there are no interactions. To define it, 
we again introduce the Green's function Gij(fl). For a 
translationally invariant system, it is a function of mo- 
mentum k. Otherwise, it has the same matrix structure 
as the noninteracting Hamiltonian Eq. (76). 

We construct the following matrix out of the Green's 
functions (here, as before, oj is taken to be real) 

Q(uj,k) = G _1 (tt«;,k)EG(tt«;,k). (84) 

Here matrix multiplication is implied, while G _1 is the 
inverse of the matrix G(Cl, k). It is clear that 



Q 2 



i. 



(85) 



We now construct the expression 



c a 1 ,a 2 ,...,a D 



/•oo f- D 

try duj J d d kQ Y[ d an Q. (86) 



Here D = d + 1 is the dimension of space-time, di = dki 
and d implies <9 W . This expression is a topological 
invariant, which mathematically reflects the fact that 
tt d (U(2M)/U(M) x U{M)) = Z for even D. One can 
check that for the noninteracting systems where G = 
[CI -U]' 1 , Eq. (86) reduces to Eq. (77). 

Note that the integration over to goes from to infinity, 
not from — oo to oo. The reason for that is that the 
integral from — oo to oo is zero, as can be verified using 



OH = EtQ(-w)E, 



(87) 



which in turn follows from the definition Eq. (84) as well 
as the chiral symmetry Eq. (12). 

To understand why Eq. (86) is a topological invariant, 
we need to further elucidate the structure of the chiral 
Green's function. Taking advantage of the explicit form 
of the matrix S given in Eq. (74), we arrive at the fol- 
lowing structure 



i=i 



(88) 



a/3\ 
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Here the matrices A, W satisfy, thanks to Eq. (12) 

A< 1,2 >(fi) = -A^' 2 \-n), W(-Sl) = W(Q). (89) 

For completeness we give the explicit expressions for 
these matrices in the appendix B. 

For a translationally invariant system, the Green's 
function can be written in the momentum space as 



G(fi,k) 



A&(n,k) w ab (n,k) 
Ai 2 6 } (n,k), 



(90) 



In the absence of any excitations around zero energy, that 
is, in case when one deals with an insulator, the functions 
A, W, are regular at O = 0. In particular this means that 
A (SI = 0) must vanish, while W(0) is a constant which 
for now we assume to be nonvanishing. This leads to 



Q(w -> 0) 



-1 
1 



(91) 



At the same time, at SI — > oo, the Green's function must 
behave as I/O, so W vanishes in this limit and 



Q(uj — > oo) 







(92) 



Thus the topological invariant measures how Q interpo- 
lates between these two values. 

Finally, one can check that Eq. (86) is indeed a topo- 
logical invariant by checking that if Q is varied slightly, 
to Q + 5Q, Eq. (86) does not change. This turns out to 
be a direct consequence of Q 2 = 1. 

The topological invariant Eq. (86) is a good starting 
point to analyze the chiral topological insulators in 3D, 
including exotic singlet time-reversal invariant supercon- 
ductors (class CI), Helium III-B (class Dili) and a parti- 
cle hopping on a bipartite lattice (class AIII) 24-26 . How- 
ever, in what follows we will use it to analyze chiral prob- 
lems in one spacial dimension. 



C. Low energy excitations at the boundaries of 
interacting one dimensional class AIII topological 
insulators 

We now have in our possession the topological invari- 
ant for interacting system in chiral symmetry classes. 
However, what we are really interested in is whether 
two insulators with two different values of I which are 
in spacial contact possess zero energy excitations at the 
boundary. More precisely, we would like to show that it 
is possible to add particles without any energy cost at 
the boundary, just like in our discussion of topological 
insulators in two spacial dimensions in Sec. V. Let us 
show that this is indeed the case. 

From now on we specialize to one spacial dimension 
(although generalizations of what follows to higher di- 
mensions should be possible to construct). Suppose we 



have a ID interacting chiral topological insulator, whose 
parameters vary in space (so that it is not translationally 
invariant). However at far away positive spacial infinity 
it becomes translationally invariant, with the topologi- 
cal invariant Ir, while at far away negative infinity its 
topological invariant is 1^. 

Just as we have already discussed in case of two spacial 
dimensions in Sec. V, as the system is no longer trans- 
lationally invariant, its Green's function is no longer a 
function of momentum only. Nevertheless, we can intro- 
duce the Wigner transformed Green's function, the one 
Fourier transformed with respect to the difference of co- 
ordinates. Such a Wigner transformed function depends 
on both the momentum k and the coordinate x. Yet at 
far away spacial infinity G becomes independent of x and 
becomes equal to its value for a translationally invariant 
system. 

We can define Q in the same way as before, 

Q{oj, k, x) = G~ (iu, k, x) S G(iu), k, x). (93) 

Here G _1 (f2,fc,a;) is defined as the inverse of the ma- 
trix G(Ct,k,x), which is not the true matrix inverse of 
G. With its help, we can define a position dependent 
topological invariant I2 by reducing Eq. (86) to the case 
D = 2 (d = 1). We denote it by I as opposed to I2 as in 
this subsection we always work in one dimensional space 
so using I without specifying dimensionality should not 
be the cause of confusion. 



I(x) 



Then 



dkQ (d u QdkQ - dkQd^Q) 



t l - i \ x =-l . J n - I \x=i 



(94) 



(95) 



for some large L such that at x — L or x — — L we are 
far from the boundary of the two insulators. 
Now we define a vector 



n a = e a i3 7 tTQdpQdjQ. 



(96) 



Here a, /3, and 7 go over values 0, 1, 2, which refer to 
the three effective coordinates, w, k and x. Notice that 

d a n a = 0. (97) 

This important relation follows from the fact that 

e«/3 7 tr daQdpQdjQ = 0. (98) 

This in turn can be proving by using that Q 2 = 1 as well 
as dQQ = —QdQ. Indeed, 

e a/ 3 7 tr daQdpQd-fQ = 

e a ^ti d a QQ 2 d /3 QQ 2 d 1 Q = 

-e a/ 3 7 tr Qd a QdpQQ 2 d 7 QQ = 

-e a/ 3 7 tr daQdpQd^Q. (99) 

We see that the expression Eq. (98) is equal to minus 
itself, so it is zero. 
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Similar manipulations can be used to prove that 
Eq. (86) is indeed a topological invariant, that is, docs 
not change under small deformations in Q, as long as 
Q 2 = l. 

Due to Eq. (97), the following integral is zero 



dS n 



0, 



(100) 



where the integral is taken over any closed surface in the 
co, k, x space. In particular, consider a surface formed by 
the planes x = L, x = —L and ui = (as well as ui = oo, 
where the integral vanishes) , with dS a being the element 
of this surface. We note that the sum of the integrals 
over x = A and x = —A is equal precisely to Ir — II- 
This allows us to express this difference in terms of the 
integral over the surface u> = 0, or more precisely 

Ir -Il= lim tr / dxdk Q (dxQ ®kQ - dkQ d x Q) . 
lbiri w->0 J 

(101) 

The limit is needed because at ui = the expression to 
be integrated is singular as we will see below. 

Now let us relate this difference to the number of zero 
energy excitations. Earlier when working in 2D we cal- 
culated the number of zero energy excitations by a clever 
application of Eq. (58). This is not possible to do in chi- 
ral systems however (formally the expression Eq. (58) is 
always zero due to Eq. (12)). Instead we take advantage 
of a similar, but distinct formula 



N = — lim ui Tr [EKd u G] 



(102) 



Here K is given by Eq. (51). We recall that as before, K 
coincides with G -1 in translationally invariant systems 
only. 

Let us check that Eq. (102) is indeed equal to the num- 
ber of zero modes. First of all, in the noninteracting case, 
G = [iui — H]~ . This allows us to write 



A = lim iuiTrY, [iui - H] 1 = Nr — N L . 



(103) 



Here Nr is the number of the right zero modes and 
Nl is the number of the left zero modes (recall that zero 

modes satisfy EV'flz, = i^jj^). This is the correct way 
of counting zero modes, since it is this difference that is 
insensitive to small changes in the Hamiltonian as dis- 
cussed after Eq. (73). 

In case when there are interactions, the analysis of 
Eq. (102) is somewhat more involved. Define the eigen- 
values and eigenvectors of the Green's function as in 
Eq. (47), or in other words 



,•(0)^(0) 



(104) 



The trace in Eq. (102) can be understood as a sum over 
these eigenvectors 



N = — lim ui ^ EKd u Gi> w , 



(105) 



where matrix multiplication is implied. 

If ^(f2) is an eigenvector of the Green's function G(f2) 
with the eigenvalue A(f2), then 'S'tp(-fl) is an eigenvector 
of the Green's function G(f2) with an eigenvalue — A(— fi), 
as a consequence of Eq. (12). 

A special role is played by the eigenvectors ifi(0). In- 
deed, if it is an eigenvector of G(0), then S^(0) is also 
its eigenvector with the eigenvalue — A(0). These eigen- 
vectors do not contribute to Eq. (105) in the limit of 
ui -> 0. 

However, there are also the eigenvalues such that 
A(fi) ~ CI as Q -> 0, or such that A(Q) ~ as Q -> 0. 
Then ^(0) and TJi/j(0) are both eigenfunctions with the 
same eigenvalue. In that case, just like in the analysis in 
the previous section Eq. (73), tp(Q) are the eigenvectors 
not only of G but also of S, with the eigenvalues +1 or 
-1. 

We call the eigenvectors with the eigenvalues A(O) ~ 
l/fi right and left poles (right poles satisfy T,ipR = tpR, 
left poles satisfy Xt/>l = —iPl), while we call the eigenvec- 
tors with the eigenvalues A(f2) <~ fi right and left zeroes. 
Then by direct substitution we find, from Eq. (105), 



N = Nr-N l -Z r + Z l , 



(106) 



where Zr and Zl are the number of right and left zeroes 
of the Green's function. 

Now we can relate A^ to the difference between the 
values of the topological invariant by following the ap- 
proach described in Sec. V. This involves expressing the 
product of two Green's functions in Eq. (102) in terms of 
their Wigner transform using gradient expansion, as well 
as expressing K in terms of G _1 using Eq. (63). This 
produces two terms. The first one, where we use the 
gradient expansion while substituting G _1 for K, looks 
like 



2i 



^TrS [-G^dxGG^d^G + G^dkGG^d^G] . 

(107) 

Note that this term should be zero after integration over 
k and x due to Eq. (62). Nevertheless we would like to 
keep it because it will make the final expression which 
will be integrated over k and x simpler. 

In the second term we use Eq. (63) for K to find 

-^TrS [G^dxGG^GkG^d^.G- 
G~ 1 dkGG~ 1 G x G~ 1 d u G] . (108) 

Combining these terms together gives 

A = limwtrS / dxd ^\G- 1 d x Gd k (G- 1 d u G)- 

J 47U L v ' 

G^dkGd^G^d^G)}. (109) 

At the same time, Eq. (101) can be rewritten, upon sub- 
stituting Eq. (93), as 



Ir — II = lim tr E 



dxdk 



(d x GG~ 1 dt c GG~ 1 - 
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dkGG^dxGG- 1 ) . (110) 

The expressions Eqs. (109) and (110) are not identical. 
Nevertheless, they are related in the following way. The 
entries of the Green's function, as introduced in Eq. (88) , 
satisfy the relations given by Eq. (89). As a result, the 
Green's function has to have one of the following possible 
asymptotic behaviors as w is taken to zero. The Green's 
function might have a pole as ui is taken to zero at this 
value of k and x leading to 

° ( "-*' x, ~(^r, <m) 

Alternatively, it is regular 

G( Wl *,*)~f^(*.*) W $f\). (112) 

Here 5^ and 8^ are the appropriate expansion coeffi- 
cients of AW and A< 2 ) from Eq. (88). 

An interesting feature of both of these expressions is 
that, in the limit to — > 0, Q is either S or — S, just as 
in Eqs. (91) and (92) and is k and x independent. As a 
result, naive substitution of either Eq. (Ill) or Eq. (112) 
into Eq. (110) gives zero in the limit of w — > 0. However, 
at the values of k and x where 5 vanish in Eq. (Ill) or 
W vanish in Eq. (112), Q varies rapidly, giving rise to 
a delta-function-like contribution to Eq. (110). Notice 
that the values of k and x where W vanishes corresponds 
to the zero of the Green's function. We mostly elabo- 
rate on this argument in two and one dimensional space, 
although it should be generalizable to higher dimensions. 

One can check by a direct substitution that if Eq. (Ill) 
holds, then Eq. (109) is equal to Eq. (110). If, how- 
ever, Eq. (112) holds then Eq. (109) is equal to minus 
Eq. (110). This was checked for G being a two by two 
matrix and a four by four matrix and we conjecture (al- 
though this has not been checked) that this is valid for 
matrices of all sizes. 

As a result, we find that 

Ir-Il=P + Z, N = P-Z, (113) 

where P is the contribution of the poles and Z is the 
contribution of the zeroes. Since the change in the topo- 
logical invariant must be an integer, P and Z are also 
integers. 

Therefore, we conclude that the change in the topo- 
logical invariant is due to either the poles of the Wigncr 
transformed Green's functions or the zeroes of the 
Wigner transformed Green's functions, which occur at 
some values of x and k, where the appropriate values of 
x are close to the domain wall separating the two topo- 
logical insulators. The poles' contributions are exactly 
opposite to those of the zeroes. 

In particular, the topological invariant will change even 
if there is no poles, P — 0, as long as there are zeroes 



present, Z / 0. Then two topological insulators with 
distinct topological invariants will have no zero energy 
excitations at their boundary. It is tempting to conclude 
that this scenario is at work for the one dimensional BDI 
systems discussed in recent papers 10-12 . However, it is 
not currently clear why the topological invariant consid- 
ered there had to change by 4 (by 8 if one works with 
Majorana chains instead of Dirac chains considered here) 
and not by 1 as would have been natural here (see ap- 
pendix C to see that Z can be any integer). 



VII. CONCLUSIONS 

We have presented a scenario in this paper where an 
interacting topological insulator can lose its edge states 
as the interactions are turned on. The method developed 
here has not yet been developed to a point where it al- 
lows to understand when this actually happens. Instead, 
it explains the mechanism by which this can happen. One 
should also be careful in interpreting the results of this 
formalism. The edge states are understood here as poles 
in the single particle Green's functions. If there is trans- 
port along the edge in such a way that the elementary 
particle carrying the current along the edge are collec- 
tive excitations, this formalism may interpret it as the 
absence of the edge states. 

In this work we restricted our attention to the topo- 
logical insulators with the topological invariant Z. The 
relationship between the invariant and the edge states 
was explored in two and one spacial dimensions only. It 
is straightforward to generalize this to higher number of 
dimensions, in particular to three dimensional space. It is 
also important to generalize this formalism to the invari- 
ant Z 2 , especially in view of current interest to 2D and 
3D All topological insulators. The expression for the Z 2 
topological invariant in terms of the Green's functions 
was found in Ref. [27]. What remains is to find what it 
tells us about the edge states (and zeroes of the Green's 
functions). All of this will be the subject of future work. 
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Appendix A: Alternative generalization of chiral 
symmetry in the presence of interactions 

Instead of Eq. (18), we could use the following con- 
struction to generalize the chiral symmetry to the inter- 



15 



acting case. Consider an operator £ such that 

&H± = -H. (Al) 

This operator does not mix creation and annihilation op- 
erators, so that its action on those operators is simply 



E^E = £,-,-a,-, E f ajE = a'-E 



*.tst.. 

3 3 % 



(A2) 



We can check that nonintcracting chiral Hamiltonians au- 
tomatically satisfy Eq. (Al) due to Eq. (70). An example 
of an interacting Hamiltonian in ID invariant under this 
operation would be (here we have spinful fermions hop- 
ping on a ID lattice) 



H = E tl + 4,i+i a ™) + ( A3 ) 

U E (4m4Ai«T* + 4i4< a 4- ia ti+i) • 

i 

Indeed, it changes sign under the sign change of every 
creation and annihilation operator on odd sites of the 
lattice. 

As a result of the symmetry Eq. (Al), for every state 
\n) there exist a state E |n) whose energy is exactly op- 
posite to that of |n). In particular, in addition to the 
ground state |0), the system should have a maximally 
excited state |0), with the energy opposite to that of the 
ground state. 

We can then define two Green's functions. One, de- 
noted by G(SY) is defined with respect to the ground state. 
The other, denoted G(f2), is defined with respect to the 
maximally excited state. One can now check that the 
sum of these two functions, 

G s (n) = l[G(il) + G(n)) (A4) 

satisfies the same relationship Eq. (12) which Green's 
functions for the chiral systems are supposed to satisfy. 
It follows that this is a legitimate generalization of the 
chiral symmetry to the interacting case, distinct from 
Eq. (18). It is possible to prove that in the absence of 
interactions Eq. (A4) coincides with the usual Green's 
function, but this does not have to be true in the pres- 
ence of interactions. We leave the question whether any 
interesting models of this type can be constructed and 
studied to further work. 



Et \n*) with one less particle than the ground state. As 
a result, the spectral decomposition becomes 

n lew — <°l hi W M 4 1°) , <°l 4^ f I"*) <»* I £ * 1°) 
^ij^ 1 ) - 2^ ~ h " 



+ e„ 



(Bl) 

Here \n) are the states with one extra particle compared 
to the ground state, while e„ = cj+ > 0, defined in 
Eq. (41). The second term in this expression can be 
transformed in the following way 

(0\aW\n*)(n*\±a t \0) = 



(0| StEatst \ n *) 

(n\ £* SjEE^ |0*) (0*| EE^atE \n) = 
£,, (0|a,|n) (n|4l0) 4,- 



SajE'E |0) 



(B2) 



To go from the second to the third line, Eq. (31) was 
used, while in the final line takes advantage of Eqs. (15) 
and (34). Finally, this gives 



Gy(fi) = 



(B3) 



E 



a t \n)(n\a}\0) £« (0| a, \n) (n\ a{ |0) E[ . 



+ 



il + e r . 



Note that this is compatible with the constraint Eq. (12). 

With the help of this expression, we can write the 
following explicit construction for matrices defined in 
Eq. (88). We introduce the matrix element 



U in = (0| di \n) . 



(B4) 



We also employ the notations Uan and which are 

equal to t/j„ as i belongs to the first or the second sub- 
lattice. Now Eq. (B3) gives 

20, 



A (1) - Vc/ (1) J7 (1)t a ^^ 

n 
n 

w afi = ^u&uvU 



2e„ 



ft 2 



(B5) 



Appendix C: Zeroes of the chiral Green's functions 
and the topological invariant 



Appendix B: Spectral decomposition of the chiral 
Green's functions 

The general spectral decomposition Eq. (40) valid for 
any Green's function can be further specialized to the 
case where chiral symmetry is present. In the presence 
of chiral symmetry, any state \n) with one extra particle 
compared to the ground state acquires a partner, a state 



Suppose the chiral Green's function takes the form 

GV) = ^ w * iw§{2) ) , (CI) 

where and are real positive functions of x and 
k while W is a complex function of x and k. This is a 
2 by 2 version of the Green's function Eq. (112). Let us 
calculate the change in the topological invariant due to 
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this form of the Green's function. Substituting this into 
Eq. (110) we find 



Ir-Il 



,. f dkdx u 2 5^S^ (d p Wd x W* - d x Wd„W*) 

hm / — ?j — . 

"^°J 2tti ( < 5(i)«5(2)a; 2 + WW*) 

(C2) 

Only the terms with nonzero ui — > limit are shown. 

Suppose W vanishes at some special values of k and x. 
At this point, in the limit oj — > 0, the Green's function 
obviously vanishes, or as we say, has a zero. 

We write W = X + iY where X and Y are real and 
they both vanish at those same values of k and x. Then 



lim 



dkdx uj 2 5^5^ {d x Xd k Y - d k Xd x Y) 



/dko 
IT 



(5W5Wuj 2 + X 2 + Y 2 y 



(C3) 



As can be checked in a straightforward manner, 



lim 



>o (cj 2 + X 2 + Y 2 )' 



= ttS{X)S{Y). (C4) 



This gives 



I R - I L = ± / dXdY6{X)6{Y) = ±1. 



(C5) 



One can check that the result of substitution of the form 
Eq. (CI) into Eq. (109) gives exactly the same expression 
as Eq. (C2) except with the opposite sign. 



1 C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

2 B. A. Bernevig and S.-C. Zhang, Phys. Rev. Lett. 96, 
106802 (2006). 

3 M. Konig, S. Wiedmann, C. Briine, A. Roth, H. Buhmann, 
L. W. Molenkamp, X.-L. Qi, and Shou-Cheng, Science 
318, 766 (2007). 

4 D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, 
and M. Z. Hasan, Nature 452, 970 (2008). 

5 S. Ryu, A. P. Schnyder, A. Furusaki, and A. W. W. Lud- 
wig, New J. of Phys. 12, 065010 (2010). 

6 A. Kitaev (2009), arxiv:0901.2686. 

7 X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 78, 
195424 (2008). 

8 G. E. Volovik, The Universe in a Helium Droplet (Oxford 
University Press, Oxford, 2003). 

9 S. Ryu, J. E. Moore, and A. W. W. Ludwig (2010), 
arxiv: 1010.0936. 

10 L. Fidkowski and A. Kitaev, Phys. Rev. B 81, 134509 
(2010). 

11 L. Fidkowski and A. Kitaev (2010), arxiv: 1008.4138. 

12 A. M. Turner, F. Pollmann, and E. Berg (2010), 
arxiv: 1008.4346. 

13 A. Altland and M. Zirnbauer, Phys. Rev. B 55, 1142 
(1997). 

14 M. Zirnbauer, J. Math. Phys. 37, 4986 (1996). 

15 M. S. Foster and A. W. W. Ludwig, Phys. Rev. B 73, 



155104 (2006). 

A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshin- 

skii, Methods of quantum field theory in statistical physics 

(Dover Publications, New York, 1975). 

A. Kamenev, in Nanophysics : coherence and transport: 

Les Houches 2004 session LXXXI (Elsevier, Amsterdam, 

Boston, 2005), arXiv: cond-mat/0412296. 

L. D. Landau and E. M. Lifshitz, Statistical Physics 

(Butterworth-Heinemann, Oxford, UK, 1980). 

D. J. Thouless, M. Kohmoto, M. P. Nightingale, and 

M. den Nijs, Phys. Rev. Lett. 49, 405 (1982). 

Q. Niu, D. J. Thouless, and Y. S. Wu, Phys. Rev. B 31, 

3372 (1985). 

A Wigner transform of a function of two variables f(x, x') 
is defined as f(k,R) = / dr e ikr f{R + r/2, R - r/2). 
V. Gurarie and L. Radzihovsky, Phys Rev B 75, 212509 
(2007). 

W. Su, J. Schrieffer, and A. Heeger, Phys. Rev. Lett. 22, 
2099 (1980). 

A. P. Schnyder, S. Ryu, and A. W. W. Ludwig, Phys. Rev. 
Lett. 102, 196804 (2009). 

P. Hosur, S. Ryu, and A. Vishwanath, Phys. Rev. B. 81, 
045120 (2010). 

M. A. Silaev and G. E. Volovik (2010), arxiv : 1005 . 4672. 
Z. Wang, X.-L. Qi, and S.-C. Zhang (2010), 
arxiv: 1004.4229. 



